Modeling quark-hadron duality in polarization observables 
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We apply a model for the study of quark-hadron duality in inclusive electron 
scattering to the calculation of spin observables. The model is based on solving the 
Dirac equation numerically for a scalar confining linear potential and a vector color 
Coulomb potential. We qualitatively reproduce the features of quark-hadron duality 
for all potentials considered, and discuss the onset of scaling and duality for the 
responses, spin structure functions, and polarization asymmetries. Duality may be 
applied to gain access to kinematic regions which are hard to access in deep inelastic 
scattering, namely for xsj 1, and we discuss which observables are most suitable 
for this application of duality. 

PACS numbers: 12.40.Nn, 12.39.Ki, 13.60.Hb 

I. INTRODUCTION 



Quark-hadron duality is a fascinating phenomenon that was first observed by Bloom and 
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and theorist _|1 

3611 communities, due to an interest in 



duality itself, and due to the huge field of experimental applications of duality in kinematic 
regions that are very difficult to access without it. 

The most straightforward definition of quark-hadron duality says that any hadronic pro- 
cess can be described in terms of either a quark and gluon picture, or in terms of a purely 
hadronic picture, provided either calculation contains all Fock states. However, in the for- 
mer case, a full numerical solution of QCD is prohibitive in most situations, despite the 
impressive progress of lattice QCD, and in the latter case, a full hadronic solution, e.g. em- 
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ploying an effective field theory, also is not feasible unless the kinematic region is restricted 
to low energies and momenta. Thus, this most general version of duality is not very useful, 
as many interesting processes take place in a region that is neither perturbative nor very 
low energy. 

There is a much more practical version of duality: in certain kinematic regimes, properly 
averaged hadronic observables can be described by a perturbative QCD (pQCD) calculation. 
This version of duality is highly relevant as perturbative QCD calculations can be performed. 
Using duality, these pQCD calculations can then be related to averaged data taken in the 
resonance region. Quark-hadron duality has been observed experimentally in many pro- 
cesses: it was discovered by Bloom and Gilman in inclusive, inelastic electron scattering, it 
made its way into the textbooks in e"'"e~ — >■ hadrons, was studied in the semileptonic decays 



of heavy mesons 



37 



38|,|39|, is considered in the analysis of heavy ion reactions '4(J], and 



forms the basis for using QCD sum rules |4l|. In addition to the "classical" examples and 
applications of duality, duality ideas are applied in new areas, too. For neutrino scattering, 
the beam energies are not well known, and an averaging will thus take place almost automat- 
ically. The application of duality is discussed for several planned neutrino experiments, see 



e.g. '4^, and duality ideas have been applied in 



43l | to nucleon/nuclear duality in neutrino 



scattering. There is also interest in dualit y in parity violation experiments 4^, and with 



regard to generalized parton distributions 



13,k5!|. A very local version of duality - assum- 

nn 



ing that it holds for just one resonance - has been used in [2a, |22| to extract information 
on structure functions at xbj — ^ 1 in the scalin g lim it from form factor data. These ideas 
were also applied to neutrino-nucleon scattering j22|. Duality ideas might also be useful for 
pion photoproduction Duality is a major point in the 12 GeV upgrade of CEBAF at 
Jefferson Lab [l^ . 

In this paper, we investigate duality in inclusive, inelastic electron scattering. New exper- 
imental data from Jefferson Lab and DESY have impressively confirmed that quark-hadron 
duality is valid down to rather low four-momentum transfers, and for many observables: 
duality in F2 was confirmed to hold down to ~ 0.5 GeV^ P|, and very recently, the 
longitudinal structure function Fl and the purely transverse Fi were separated, and found 
to exhibit duality for > IGeV^ j3|. Experimental evidence for duality in spin observables 
has been reported for from Hermes 0] and from Jefferson Lab for the first moment 
of Qi- While this is exciting all by itself, these data have inspired experimentalists to apply 
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duality to the extraction of information on the deep inelastic region in kinematics that are 
not readily accessible. Duality allows us to connect the perturbative regime of quarks and 
gluons with the strongly non-perturbative resonance regime. The earliest example discussed 
was the extraction of the elastic nucleon form factor from the deep inelastic scaling curve 

n n 

j47| . In |4Sl| , higher twist contributions were inferred from the resonance data. 

These two kinematical regions have traditionally been separated, as it was believed that 
the physics of quarks and gluons had little connection to the collective phenomena of reso- 
nances. At low invariant masses W of the final state in an (e, e') reaction on the nucleon, 
one observes many resonance bumps in the cross section, and W < 2 GeV is traditionally 
referred to as the resonance region. For higher invariant masses, W > 2 GeV, the cross 
section becomes smooth and exhibits Bjorken scaling, and this region is referred to as the 
deep inelastic region. Note that this strict division of kinematics was introduced historically, 
even though non-resonant processes are present for low W, and resonances with larger mass 
may contribute for high W. The demarcation line of W = 2 GeV is plotted in the xbj - 
plane in Fig. 

One consequence of this traditional subdivision is that large amounts of data with W < 
2 GeV were cut from deep inelastic analyses of data, leading to a paucity of data at very high 
Bjorken xbj- The region of large x^j, xbj — ^ 1, is referred to as the (deep) valence region, and 
is the subject of much interest. In particular, one would like to study the valence quark spin 
distribution of the nucleon. This can be achieved by measuring the polarization asymmetry 
Ai of the proton and neutron. For xbj 1, many, and widely different predictions, exist for 
the polarization asymmetry of the neutron, running the gamut from 0, predicted in unbroken 



SU(6), to 1, predicted in pQCD, and everything in between, see j49|,|50| for good reviews of 
the situation. In the chiral soliton model, there are even predictions of negative values Isi]. 

The experimental data available for at high xbj was rather scarce and afflicted with 
very large error bars for xbj > 0.4 before the advent of the recent Jefferson Lab data 
They also extended only up to xbj ~ 0.6 [3, 1^ 1^ H]- This is due not only to the 
fact that polarization experiments are always more difflcult to perform than unpolarized 
measurements, but mainly due to the fact that in order to achieve xbj —>■ 1 for the deep 
inelastic regime, one has to use very high four- momentum transfers Q^, see Fig. ^ This 
drastically reduces the cross section, as the cross section is proportional to the Mott cross 
section, which in turn is proportional to Thus, accessing one and the same xbj in 
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FIG. 1: Plot of the XBj and kinematic plane. This kinematic plot shows the deep inelastic 
regime and the resonance regime. The solid line corresponds to an invariant mass of the final state 
of = 2 GeV. Every point above the line lies in the resonance region, W < 2 GeV, and every 
point below the line lies in the deep inelastic region W > 2 GeV. The invariant mass of the final 
state, W, is related to the four momentum transfer and Bjorken xbj by W"^ = + Q^{^^ — 1), 
where M is the nucleon mass. 

the resonance region and in the deep inelastic regime leads to much lower count rates in 
the deep inelastic regime. For example, a measurement at xbj = 0.8 can be performed at 
Q"^ = 2 GeV^ in the resonance region, or at > 15 GeV^ in the deep inelastic region. The 
count rate in the deep inelastic region will be lowered by more than a factor of 50 compared 
to the resonance region measurement. This makes taking data in the deep inelastic regime 
for large xbj extremely difficult. A recent Jefferson Lab experiment improved the situation 
by measuring in the deep inelastic regime, up to xbj > 0.6 with very reasonable error bars. 



50,^ 



thus decreasing the uncertainty by an order of magnitude compared to older data 
Still, the experimental exploration of at really large x^j — > 1 has not been feasible yet, 
thus making it the prime application for duality. Even with the planned 12 GeV upgrade 
of Jefferson Lab, data in the deep inelastic regime will be accessible only up to xbj ~ 0.75 

Q. 

As the kinematics plot shows, large values of xbj can be accessed at low in the 
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resonance regime. A measurement of - or any other observable of interest - in the 
resonance region can then be averaged, and will yield the same information as a direct 
measurement in the deep inelastic region, provided that duality holds. Using this approach, 

could be measured up to xbj ~ 0.9, by taking data in the Delta resonance region j^. 

Another interesting application of duality that was recently discussed is the application 
of duality to the EMC effect j^. Even though the EMC effect has been mapped out for a 
large kinematic region, data at very high xbj and for lighter nuclei are scarce. The data 
base could be significantly expanded by the application of duality to lepton scattering from 
nuclei. 

Before we can embark on this new experimental approach, we need a good, solid under- 
standing of how well duality holds, and where it holds. Ideally, one would come up with a 
certain accuracy of the duality procedure, which then could be quoted as a systematic error 
for the extraction of A" or other observables in the relevant region. Here, we need more 
theoretical input and guidance. As we currently do not understand quark-hadron duality 
from first principles, modeling is our best tool to obtain the answers necessary for applying 
duality to extract A". 

Currently, theorists tackle this problem by modeling duality in two different ways: one 
branch starts out from the non-relativistic constituent quark model, with some relativistic 



17, 



1£ L the other branch starts the modeling 
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The former branch 



corrections, to describe duality [11|, ll^ 
with a relativistic one-body equation 
makes contact with the phenomenology. It was started by the pioneering work of Close 
and Isgur h|, where the authors investigated how a summation over the appropriate sets 



of nucleon resonances leads to parton model results for the structure function ratios in the 



12 1 to include the effects 



SU(6) symmetric quark model. This work was recently expanded 
of SU(6) spin- flavor symmetry breaking. In the authors considered the first five 

low-lying resonances, and performed a careful analysis of the onset of duality for F2 and 
Qi. Our results belong to the latter branch. The goal of these modeling efforts is obvious: 
to gain an understanding of quark-hadron duality and the conditions under which it holds, 
by capturing just the essential physical conditions of this rather complex phenomenon. We 
imposed these basic requirements for a model: we require a relativistic description of confined 
valence quarks, and we treat the hadrons in the infinitely narrow resonance approximation. 
This paper is the fourth in a series of papers, in which we have modeled duality with 
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increasing complexity. All models that we have presented so far have reproduced the features 
of duality in a qualitative manner. We started out with an all scalar model and 
gradually improved the model until all the particles had proper spin we focused 

on the model results for the unpolarized responses. For the first time, we investigated the 
dependence of our results on the type of potential we employed. In the present paper, we 
focus on the spin observables: the responses Rt' and Rtv that are accessible only with spin, 
the spin structure functions gi and (72? and the polarization asymmetries Ai and A2. 

This paper is organized as follows: in Section |nl we briefly state the properties of the 
model, then, in Section ITTH we discuss some of the formalism for polarized inclusive electron 
scattering. In Section IIVI we discuss our numerical results for the various spin observables 
for the bound-free and the bound-bound transitions. We end with a brief summary of our 
results and an outlook. 



II. THE MODEL 

We use the same model as in js^. For the convenience of the reader, we present the key 
ingredients of our model here. 

Our model consists of a constituent quark bound to an infinitely heavy di-quark and is 
represented by the Dirac hamiltonian 

= a-p + /?(m + K(r)) + K(r), (1) 

where the scalar potential is a linear confining potential given by 

K(r) = 6r, 6 = 0.18GeV^ (2) 

We have used the constituent quark mass in this paper, as our main interest is the study 
of quark-hadron duality, which sets in at rather low Q^, experimentally ~ 0.5 GeV"^ 
is enough. In this kinematic region, the appropriate degree of freedom is the constituent 
quark, which has acquired mass through spontaneous chiral symmetry breaking. We have 
used a value for the quark mass of m = 258.46 MeV - obtained previously in a fit to heavy 
mesons jssj. However, nothing hinges on using that particular value: we changed our quark 
mass to m = 10 MeV , in order to have a value reminiscent of a current quark mass, and 
repeated our calculations. It turns out that, while scaling does set in a little faster, there 
are no qualitative changes in the results. 



7 



In our model, the vector potential is provided by a vector color Coulomb potential. 
Calculations will be presented where the vector color Coulomb potential is absent, that is 
Vv{r) = 0, where the vector potential is the simple static Coulomb potential 

= K(r) = -^^ (3) 

with as = 0.181 and where the color Coulomb potential is corrected to allow for the running 
coupling constant in a manner similar to that used by Godfrey and Isgur [sol. The vector 
potential then has the form 

K(r) = VUr) = -- ( + 5^a.erf(7.r) j (4) 

where 





= 0.118 


Po 


= 0.04 


6 


= 0.01 


ai 


= 0.239 


a2 


= 0.271 


7i 


= 0.746 GeV 


72 


= 5.40 GeV. 



(5) 

Note that we use different scalar and vector potentials, in contrast to |21|,|2^, where K = K 
is used to simplify the calculations. 

We assume that only the light quark carries a charge, and we choose unit charge for the 
light quark for simplicity. 

III. SPIN OBSERVABLES 

In this section, we briefly review the formalism for calculating responses for targets with 
arbitrary polarization axes, and connect the definitions of the polarization asymmetries 
Ai,A2 and the spin structure functions gi,g2 to the responses. 

The hadronic tensor for targets with an arbitrary polarization axis in s direction is: 

= E <m" < + S-s)\\m> (6) 
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with 

uj''^^„= <lQ]-m\r^\n'l'j'm' ><n'l'j'm'\r\l{}]-m'' > (7) 

n'V j'm' 

where J is the electromagnetic current operator. The ground state's z-component of j is 
denoted m instead of nij for brevity. Components of the hadronic tensor can be combined 
to give 



Wl 


= lyoo 


Wt 


= + W 


Wtt 


= 2^{W+^) 


Wtl 


= -23?(H^°+ 


Wtv 


= -23?(W^°+ 


Wt' 





(8) 

The cross section for electron helicity h and target polarization axis s can be expressed in 
terms of the hadronic tensor and the leptonic coefficients vk, K = L,T, TT, TL, T', TV as 

^^^^Q, = ctm (vlWl + vtWt + vtlWtl + vttWtt + h [vtl'Wtu + vt'Wt']) (9) 

where aMott is the Mott cross section, q is the three-momentum transfer from the electron 
to the target, v is the energy transfer and = — z/^, and h denotes the electron helicity. 
The leptonic coefficients are given by 



Q 



4 



■)2 



Q' 2 ^ 

Vt = + tan - 



2q 



, e 9 

vt' = + 2 ^^"^ 2 

For arbitrary target spin, all six combinations of the hadronic tensor are non-zero. How- 
ever, for spin 1/2 targets, only four combinations contribute: L,T,TL' and T'. Inserting 
the results for the current matrix elements, and exploiting selection rules and symmetry 
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relations between various current matrix elements, one finds for our case 



Wl 
Wt 



l/2-m^ 



m n'l'j' 



m n'l'j' 



< n'l'j'm + 1| J+|10-m > p + | < n'l'j'm - 1\J \10-m > 



/2-m\ 



m n'l'j' 



< n'l'j'm + 1| J+|10-m > 



< n'l'j'm — 1\J \10-m > 



Wtt 
Wtl 
Wtu 








2s, y m< l{)--\Jl\n'l'j'- >< n'rj'-\J+\10- --> 
X y 22' °' 2 2' ' 2 2 



(11) 



n'l'j' 



Explicitly carrying out the summation over the initial spin m for the T' and TV combina- 
tions, and substituting for the hadronic tensor combinations in the expression for the cross 
section Eq. Q, we find: 
da 



dE'dVt' 
+ h 



m n'l'j' m n'l'j' 

-2..t;TL'$^5R(J°'t(i)J+(-l))+.;^,..5^ (|J+(i)p-|J-(-i^^ ^12) 

n'l'j' n'l'j' ^ ^ J ^ 

where we abbreviated the current matrix elements as J^{m) for < n'V j'm\J^\lQ}^m >. 
Now we can relate this cross section to the definition of the polarization observables. The 
polarization asymmetry A\\ is defined as the ratio of the difference and sum of the cross 
sections for longitudinally polarized electrons and target polarization parallel or anti-parallel 



to the beam, see e.g. 



laiiy 

SI 



da- 



da^ 



dE'dQ' dE'dQ' 



I 



da-^ da^ 
+ 



dE'dn' dE'dQ' 



(13) 



The upper superscript denotes the direction of the electron polarization, the lower super- 
script indicates the direction of the target polarization. The relevant cross sections are 
obtained by using /i = 1 in each case and s and —s for the target polarizations parallel and 
antiparallel to the beam. Using the conventional coordinate system, this means that the tar- 



get is polarized in the x — z plane, see e.g. 



61 



In this case, s = (sin a, 0, cos a), where a 
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is the angle between transferred momentum q and beam momentum k: cos a = +^^beam'^ _ 
Substituting and rearranging now yield: 

Now, we compare this expression to the definition of the polarization asynimetries Ai,A2. 
For the convenience of the reader, we quote the standard definitions 0,0]: 

Aii=D{A^+r]A2) (15) 

with the depolarization factor D 

where y = and i? = ^ = ^(1 + ^) - 1, and 

with the magnitude of the virutal photon's longitudinal polarization, e = (1 + 2^ tan^ ^)~^ 
and 7 = ^^^f£j . Note that the ?/ defined here is not the variable used in the next section. 
Thus, we can read off the expressions for the polarization asymmetries in terms of response 
functions from Eq. (fT^ as 



R 



where we used the symmetry of the current matrix elements and the definitions for the 



responses 



7-0 12 



R^^, = _2^5fJ(jO't(j+ + J-)) (19) 
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The symbol indicates the average over initial states and the sum over final states. The 

spin structure functions can be found as functions of the responses using the relation between 
polarization asymmetries and spin structure functions: 

A. = 7^ (20) 
with the unpolarized structure function Fi = MWi = ^MRt- We find: 

g, = -1m4(/?t' + ^— /2tl') 
2 V2 qi^ 

92 = l^^^iRT' - -^^Rtl') (21) 

Summarizing our expressions for the polarization observables in terms of response functions, 
we have 

Rt' 



^1 
A2 



Rx 

1 QRtl' 
V2q Rt 



g, = --M-{Rt' + -j=^Rtu) 
2 q^ ^qv 

g^ = ^M'^iRT. - ^-Rtu) . (22) 
2 q^ V2 g 

For completeness, we also quote the expressions for the unpolarized structure functions 
Wi and W2 in terms of the responses: 

W,{Q\u) = ]^RT{q,u) (23) 



and 

W2{Q\ v) = ^RL{q, y) + ^Mq. y) • (24) 



-^RL{q,-) + ^2- 



We focused on the longitudinal and transverse response functions in our last paper js^, 
and we now focus on the true spin observables, which we can now calculate, as quark spin 
is included in our present model. 

The Dirac wave functions and energy eigenvalues are obtained by integrating the Dirac 
equation using the Runge-Kutta-Feldberg technique and solutions are obtained for energies 
up to 12 GeV with the radial quantum number of n = 200 and \k,\ < 70. 
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ow it to decay. We refer to this process as the bound-bound transition 
3ol |. we referred ot it as FSI). Thus, we do not include any particle production 



In our model, we excite the bound quark from the ground state to higher energy states, 
and do not al 
(previously, in 

in our model, and are strictly quantum-mechanical in this sense. We do not have any gluons 
in our model, either, which means that we do not encounter any radiative corrections. Since 
the response functions consist of a sum of delta functions, we choose to smear out the 
response functions by folding with a narrow gaussian for purposes of visualization. The 
smeared response functions are then given by 



oo 

RK{q, y) = ^ I rfz/'e-'^i?r™'(g, , (25) 
V^e J-oo 

where K stands for L, T, T' or TV . 

Before presenting numerical results, we would like to remind our readers that, while the 
present model is more realistic than its predecessors, its results should not be compared 
quantitatively to inclusive electron scattering from a nucleon. Due to the assumption of an 
infinitely heavy antiquark (or diquark) to which the light quark is bound, our calculation 
most resembles inclusive electron scattering from a B-meson, which has never been measured. 
The goal of our work is to gain a qualitative understanding of duality, and the current 
simplification is no impediment to this. 

As discussed in ^|, the current matrix elements naturally are functions of the three- 
momentum transfer q. Therefore, it is convenient to show our results for fixed three- 
momentum transfer g as a function of the y-scaling variable 



y 



^{v + Eo)2 - m2 - q . (26) 



The physics are not affected by presenting our results in this fashion. In fact, later on we 
will exphcitly compare our results plotted for fixed q and y to our results plotted for fixed 
four-momentum transfer and the appropriate x scaling variable. 

In the following, we will present analytic and numerical results for the bound-free and 
bound-bound transitions, and we investigate if duality holds or not. The conditions that 
need to be fulfilled to see duality are: scaling of the bound-bound transition and bound- 
free transition to the same scaling curve, and oscillation of the bound-bound results at low 
momentum transfer around the scaling curve. We are going to check if our model results 
qualitatively reproduce the signature of duality as seen in the electron scattering data. 
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In each case, we will investigate the approach to scaling for the different observables, the 
scaling curves themselves, and the behavior at low momentum transfers. 



Now, we turn to the numerical results of our model calculations. Within our model, we 
calculate two different processes: the bound-bound transition and the bound-free transition 

n 

of the light quark. The bound- free transition (referred to as PWIA in |30^]), is the analog 
to perturbative QCD. The bound quark is knocked into the continuum by the absorption of 
the virtual photon. 

We discuss the bound-free transition first, employing just the linear confining potential. 
We investigate the scaling results and the onset of scaling for the bound-free transition. 
Next, in Section riV(^[ we show our results for the bound-bound transition and the linear 
confining potential. There, we focus on the onset of scaling and the low-g duality. Then, 
we take a look at the role of the ground state p-wave contribution in Section IIVDI Finally, 
we investigate the effect of employing a static Coulomb and a running Coulomb potential 



For the bound-free transition, we can reach arbitrarily high energy transfers without any 
numerical complications, as we do not have to solve for any high energy bound states in the 
final state. Moreover, we can determine the analytic expressions for the responses, and thus 
for all other observables, for g — > oo at fixed y. As the responses Rl and Rt-, which can be 
accessed in unpolarized scattering, enter some spin observables, we quote all four responses 
for the reader's convenience 13011: 




IV. NUMERICAL RESULTS 



in Section HV El 



A. y-scaling for the bound-free transition 




167r2, 



— — / dppl^{y + qY + ni^ril{p)+mns{p) 




(27) 
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1 /•!/+2g r 

RTiq,y) = / dpp<^y{y + q)^ + m'^n^^{p)-mns{p) 

1 J\y\ I 



+ 2qy + - + 2qy - p 



2q 



2pq 



Kip) > , 



(28) 



RT'{q,y) 



8n^q 



y+2q 



dpp< [\/{y + q)'^ + m'^ + m] n_{p) - pnl{p) 



2\ 2 



2pg 



+ + + ^s(p) -mnl{p) - ^ — —Kip) } > 



and 



RTL'iq,y) 



J/+2? 



2p 



^ _ 1^ (y^ + 2yq - p^ 
2pq 



(29) 



n-ip) + qrisip) 



y^ + 2|/g - 
2pg 



{mnl(p) -pusip)) 



In the limit of large g, the bound-free response functions become 



lim Riiq.y) 



y 

. dpp'tnlip) + -KiP) 



(30) 



(31) 



and 



lim Rxiq^y) = ^ dpPS KiP) - - < 



lim RT'iq,y) = „ 2 



1 r , /</' 



n_(p) +ns(p) - -<(p) } , 
P 



(32) 



(33) 



\/2 r°° ( v"^ I 1 

lim i?TL'(g, y) = ^ dpp\{l-^)- n_{p) + n,(p) ^ . (34) 
These response functions therefore scale in y. 

The vector and scalar momentum density distributions n„(p) and Usip) are defined in 
terms of the ground state wave function 

>'!oj(p)^ri("p) 

as 



, , , , rziip)y7d^v) 

^2 ^2 



(35) 



^j;(p) = (nlip), A<(p) 



PI 



(36) 
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with 
and 

<(p) = ^^;;l(p)4i(p); (38) 

and 

ns{p) = ^{^^^,r{p)-^[-f{p)) ; (39) 

moreover, 

n_(p) = l(^;;j'(p)) . (40) 

The asymptotic responses for the bound-free transition are shown in Fig. |2l The two 
purely transverse responses, Rt and Rt', have opposite signs but similar peak positions at 
slightly negative y values and peak heights. The longitudinal response i?^ is the smallest 
of the four responses, and the only one to peak at a slightly positive y. The transverse- 
longitudinal interference response is the largest response, and peaks aX y = 0. This response 
is symmetric around y = as it contains only terms proportional to y^ and y'^, but no 
contributions linear in y. 
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FIG. 2: The asymptotic behavior for g ^ 00 and fixed y of the responses for the bound-free 
transition. We show Rj^ (solid line), Rt (dotted line), Rti (dashed line), and Rtl' (dash-dotted 
line) . 
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The similar size and shape of Rt and Rt', together with their opposite signs, leads to a 
value of the polarization asymmetry Ai = —j^ close to 1, see Fig. 01 We show Ai only for 
y values for which the responses have appreciable values, i.e. for — 1 < y < 0.6. While we 
could calculate Ai in a region where the responses are tiny, the results would be meaningless. 
Our result is close to the pQCD prediction of 1 for Ai in the limit of x — > 1. For large Q^, 
large negative values of y - much bigger than y = —IGeV - correspond to x = 1. 
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FIG. 3: The asymptotic behavior for g ^ cxd and fixed y of the spin structure functions g\ and 
52 and the polarization asymmetry A\ for the bound- free transition. We show g\ (solid line), 52 
(dotted line), and A\ (dashed line). Note that A2 vanishes in this limit. 



The asymptotic values of the spin structure functions g\ and (72 are also shown in Fig. |31 
The spin structure functions have peaks at small, negative y values. They have opposite 
signs. Of the two functions, has the slightly larger peak value. Note that the asymptotic 
value of A2 is zero. 



B. Approach to scaling in the bound-free transition 

Now, after showing the asymptotic values, we will consider the approach to scaling in the 
bound-free transition. Previously, we have studied the approach to scaling of the unpolarized 
longitudinal and transverse response functions jsO]. There, we found that the onset of scaling 
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is not influenced very much by the spin of the target particle. The two response functions 
Rt' and Rtu are accessible only with polarized beams and targets. It is interesting to check 
if they scale just like their unpolarized counterparts. 
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FIG. 4: The transverse-primed response function Rt' is plotted versus y for several low (bottom 
panel) and high (top panel) values of the three-momentum transfer q. The results shown have been 
calculated for the bound-free transition, using the linear potential. 

The scaling behavior of the transverse-primed response Rt' for the bound-free transition 
is shown in Fig. |3J The lower panel shows the response for several q values up to 10 GeV. 
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With increasing momentum transfer, the peak of the response moves to more negative values 
of y, and increases in height. The width of the response also increases. The change in the 
response when going from q = 2GeV to q = 4GeV is significant, and the difference between 
the q = 4GeV and q = 6GeV results is noticeable. Increasing q leads to very small changes, 
visible mainly at the peak and the large, negative y flank. The top panel, displaying the 
response for various higher values of q and the asymptotic value discussed above, shows that 
the response has converged to the asymptotic value roughly at g = iOGeV. 
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FIG. 5: The transverse-primed response function Rtu is plotted versus y for several values of the 
three-momentum transfer q. The results shown have been calculated for the bound-free transition, 
using the linear potential. 



The transverse-longitudinal primed response Rtl' for the bound-free transition is shown 
in Fig. The figure shows the response for lower values of the three-momentum transfer q. 
While there are small changes, amounting to a small shift of the entire response towards the 
negative y values, one can see that this response scales much faster than Rt'. We do not 
include another panel with Rtu calculated for higher q values, as the curves all coincide. 

From this, we can see that the transverse-primed response Rt' shows a scaling behavior 
very similar to Rt (as discussed in and that the only interference response that is 

accessible in inclusive electron scattering, Rtl', exhibits the earliest onset of scaling. 

Now, we will consider the y-scaling behavior of the spin structure functions. The spin 
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FIG. 6: The spin structure function g\ is plotted versus y for several low (bottom panel) and high 
(top panel) values of the three-momentum transfer q. The results shown have been calculated for 
the bound-free transition, using the linear potential. 



structure functions are simply combinations of the spin-dependent responses multiplied with 
some kinematic factors, see Section 11111 We have already seen how the relevant responses, 
R^i and -Rtl'? scale. Taking a look at the kinematic coefficients shows us that the onset 
of y-scaling will be determined by the q values at which the kinematic factors reach their 
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asymptotic values. For the convenience of the reader, we recall Eq. ()22p here: 

2 \ ^/2qv ) 2 q^ 2V2 r 

1 z/2 / \ V \ 1 z/2 1 z/3 

g2 = -M— Rt' - -^-Rtl' = -M—Rt ^M—Rtl' 

2 q^ \ ^/2q J 2 q^ 2^2 q^ 

The asymptotic values of the kinematic factors multiplying the responses are, for fixed y 
and large q, 



v'' . 1 



q^ q 



1 + -2{y - Eo) 



(f" q 



2(^0 - y) ^ 



1 + -3(y - ^o) ^ 1 (41) 



q 



leading to the asymptotic values of the spin structure functions: 



(7i ^ --MRt, 

92 ^ ]^M{Rt' - -^Rtv) ■ (42) 

From Eq. ()41|) . one can see that the relevant scale for the onset of scaling is given by 
the ground state energy Eq. All responses peak roughly in the region y ^ 0, so that the 
three-momenta q necessary to reach the asymptotic values are determined by Eq alone. The 
spin structure function gi for the bound-free transition is shown in Fig. IHl The lower panel 
shows the spin structure function for several q values up to 10 GeV. One can see clearly 
that gi changes significantly for each increase in q, and that convergence has not set in at 
q = lOGeV. This behavior is due to the slow approach of the kinematic factors in gi to 
their asymptotic values. The responses themselves scale much faster, and are very close to 
their asymptotic values at g = 10 GeV, whereas the spin structure function gi is off its 
asymptotic value by almost 50%. Results for g2 at lower (bottom panel) and higher (top 
panel) three-momentum transfers are shown in Fig. [7| Just like g2, it approaches its scaling 
limits only at high momentum transfers. 
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FIG. 7: The spin structure function g2 is plotted versus y for several low (bottom panel) and high 
(top panel) values of the three-momentum transfer q. The results shown have been calculated for 
the bound-free transition, using the linear potential. 

The polarization asymmetries Ai and A2 are ratios of response functions, and may there- 
fore show a different scaling behavior than the responses themselves. As we have seen that 
the two purely transverse responses, Rt and Rt'-, have a very similar approach to scaling, 
one may expect to see an even more rapid scaling of Ai, the ratio of the two transverse 
responses. The results for Ai are shown in the top panel of Fig. |H| As Ai is a pure ratio 
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of two responses that scale at reasonable values of g, without any kinematic factors, we can 
expect Ai to scale fast. Indeed, one clearly sees that in the region of —0.5 < y < 0, where 
the responses are largest, Ai scales almost immediately. Only the q = 2 GeV curve differs 
very slightly from all the other curves, including the asymptotic value, in this y interval. 
Outside of that region, where the responses have less strength, the scaling takes a bit longer. 
In all cases, q = lOGeV is very close to the asymptotic value. 

The approach to scaling of the polarization asymmetry A2 is shown in the bottom panel 
of Fig. IHl Note that — for g — > 00 and fixed y, because the kinematic factor ^ goes 
to zero in this limit and the responses scale in y. Thus, the approach to scaling of A2 will 
mainly be determined by the kinematic factor ^. So, even though both Rtu and Rt scale 
at reasonable values of q, the scaling of A2 is delayed. The differences between the curves 
for increasing q values decrease, but the result for q = lOGeV is still far away from the 
asymptotic value of 0. Just as for the spin structure functions, the kinematic factor that is 
present determines the very slow onset of scaling. 

In summary, we observe that the response functions reach values very close to their 
asymptotic values around g = 10 GeV, and scale at the latest at g = 40 GeV. The 
response Rtl' is the fastest scaling response, and the response Rt' is the response that 
scales most slowly. The spin structure functions and the polarization asymmetry A2 scale 
only for much higher momentum transfers, due to the kinematic factors in their definitions. 
The polarization asymmetry Ai, on the other hand, is the observable for which we observe 
scaling at the lowest momentum transfers. 

Note that the scaling behavior of the discussed observables does not change when we use 
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an "x-type" scaling variable, as in our previous papers 
Fig. ini where we show gi as a function of m = 2^(a/z^^ + — z/)(l + \J'^ + , for fixed 
Q^. Note that u g oes to = ^^^-^^j^f^xsj for large Q"^, its properties have been discussed 
previously 2^ |2^. This is the appropriate scaling variable for use with fixed even at 
low Q2. 

We have generated the plot by starting out with our results as functions of g and y, 
then calculated the corresponding values of the four-momentum transfer and the scaling 
variable u. The results were then sorted into bins. This is why there are no continuous 
lines in the figure, just single data points. 

One can see clearly that the scaling behavior of gi is independent of the variables chosen: 
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FIG. 8: The polarization asymmetries Ai (top panel) and A2 (bottom panel) are plotted versus y 
for several values of the three- momentum transfer q, as well as their asymptotic value for q ^ 00. 
The results shown have been calculated for the bound-free transition, using the linear potential. 

just as seen in Fig. IHlfor fixed q and the variable y, at low Q^, gi starts out negative, then 
changes sign, and takes a very long time to scale. Even between the points for = 16 GeV^ 
and = 20 GeV^, there is a significant difference. Thus, we see that u-scaling is just as 
slow for Qi (and g2 and A2) as y-scaling. The scaling behavior does not change at all with 
the use of different scaling variables. 
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FIG. 9: The spin structure function gi is plotted versus u for several values of the four-momentum 
transfer squared Q^. The results shown have been calculated for the bound-free transition, using 
the linear potential. 

C. Scaling and low q duality in the bound-bound transition 

Now we proceed to discuss the behavior of the responses, spin structure functions, and 
polarization asymmetries for final states consisting solely of resonances. We refer to this 
case as the "bound-bound" transition. The wave functions for the excited states have to 
be obtained by solving the Dirac equation numerically. Due to the involved nature of the 
numerics, the highest accessible momentum transfer is g = 10 GeV. Work on extending our 
calculations up to 30 GeV using a WKB approximation will be reported elsewhere. Here, 
we focus on the approach to scaling and on duality at low q. 

We start again by considering the two polarized responses. We show Rt' in the top 
panel of Fig. ^| The resonance bumps at lower q values give way to smooth curves formed 
by many, closely spaced resonances. The graph shows that there is still a small difference 
between the results for q = 8 GeV and g = 10 GeV, we have not yet reached the scaling 
value. This is to be expected, as the bound-free transition scales only at higher q values for 
this response. In general, the more complicated interplay of the various final state resonances 
leads to slower scaling for the bound-bound transition. The resonances oscillate around the 
smooth curves for higher momenta q on the positive y flank, but are a bit below the smooth 
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curves at the negative y flank. We see a very similar behavior for Rtv , in the bottom panel 

of Fig. uni 
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FIG. 10: The transverse-primed response function Rx' (top panel) and the transverse-longitudinal 
primed response function Rtl' (bottom panel) are plotted versus y for several values of the three- 
momentum transfer q. The results shown have been calculated for the bound-bound transition, 
using the linear potential. 



For Ai, the bound-bound transition results scale very quickly in the region of the response 
peaks, —0.5 < y < 0. This is consistent with the scaling behavior of Ai as observed in 
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the bound-free transition, and stems from the definition of Ai as ratio of two fast-scahng 
responses. It is interesting to note that the resonance bumps that are visible at g = 2 GeV 
have vanished almost completely at g = 4 GeV, even though they are clearly present at 
g = 4 GeV for the transverse and transverse-primed responses, whose ratio forms Ai. This 
exemplifies the very quick onset of scaling for this observable, which we can observe even in 
the region of low momentum transfer q. 

Our current model gives at best a very qualitative insight on the workings of duality for 
proton targets, as it is most closely related to electron scattering off a heavy-light meson, i.e. 
a B-meson. Nevertheless, our results seem to be very encouraging as far as the prospects of 
applying quark-hadron duality to extracting deep inelastic information on Ai from data in 
the resonance region are concerned. 
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FIG. 11: The polarization asymmetry Ai is plotted versus y for several values of the three- 
momentum transfer q. The results shown have been calculated for the bound-bound transition, 
using the linear potential. 



The polarization asymmetry A2, which contains a kinematic factor multiplying the ratio 
of two responses, is shown in Fig. ^| In contrast to Ai, A2 changes significantly with 
changing momentum transfer, and approaches its scaling value of zero only slowly. As 
mentioned when discussing the bound-free transition results for A2 in the previous section, 
the behavior of this polarization asymmetry is mainly determined by the kinematic factor. 
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not by the ratio of the responses. 

The same scenario - large variation for changes in momentum transfer q and slow approach 
to the scaling values - is repeated for the spin structure functions gi and g2- As for y42, 
their behavior is determined by the slowly-scaling kinematic factors, not by the fast-scaling 
responses. For the lowest considered momentum transfer, gi even has a different sign, 
compared to the higher momentum transfer results. 

Our results suggest that the application of duality to resonance region data for Ai should 
be quite safe, while A2, gi-, and g2 seem to be observables that are not amenable to such an 
extraction procedure. 



As far as duality in the spin structure function gi is concerned, the available data 



are not vet conclusive, and may imply a different onset of duality for g\ and g^. The 



Hermes data jS 
recent analysis 



iity 

i 



for g\ indicate that the onset of duality takes place at LSGeV^. A 

of Jefferson Lab data for the first moment of (7" indicates that duality 
could hold at as low as 1 GeV^. A dedicated Jefferson Lab experiment to study dualit 
in neutron spin structure functions in the resonance region is currently being analyzed 
Dong et al. [l^, 13, ll^ recently used a constituent quark model for the five lowest-lying 
resonances and various parametrizations of scaling data for a careful theoretical analysis of 
duality in F2 and gi for the proton target. They included target mass corrections and find 
that duality in gi is not seen below Q"^ 2GeV'^. Close and Isgur using an SU{Q)- 
symmetric, constituent quark model, also predicted a slower onset of duality for gi, and 
pointed out that duality sets in faster in this model for g^ than for g^. Dominant magnetic 
interactions are necessary for the symmetric quark model, and for the neutron target, this 
is realized. 

Our simple model is, of course, far from describing actual electron scattering measure- 
ments. It is not intended to provide a quantitative description of data, but to provide 
qualitative insights into Nature. However, our findings on the scaling behavior of the polar- 
ization asymmetries and spin structure function might provide more general guidance: what 
we see is that the onset of scaling is driven mainly by the kinematic factors multiplying the 
reasonably fast-scaling response functions in the definitions Eq. ()22|) of the spin structure 
functions and polarization asymmetries. Where the factors are missing, i.e. for Ai, con- 
vergence is rapid and the prospects for a successful application of duality are very good. 
Where the factors are present, and take a very high momentum to reach their asymptotic 
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value, e.g. in gi, scaling is very slow and duality does not hold. This suggests that the 
validity of quark-hadron duality in spin observables may be strongly related to kinematic 
factors, and less to dynamics. The faster scaling for might be explained due to a smaller 
longitudinal-transverse interference response Rtu, see Eqs. fl41f42j] . 

However, one has to be cautious as certain features of nature are not included in our 
model: it is well-known that the A resonance has an unnaturally small coefficient multiplying 
the leading term and therefore gives rise to peculiarities in the scaling and duality behavior, 
see 16^ . 
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FIG. 12: The polarization asymmetry A2 is plotted versus y for several values of the three- 
momentum transfer q. The results shown have been calculated for the bound-bound transition, 
using the linear potential. 



D. The role of the p-wave contribution 

One interesting question is which role the p-wave in the ground state plays. We will 
discuss the role of the ground state p-wave contribution for the onset of scaling for the 
observables, and for the scaling results themselves. 
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FIG. 13: The spin structure functions gi (top panel) and §2 (bottom panel) are plotted versus y 
for several values of the three-momentum transfer q. The results shown have been calculated for 
the bound-bound transition, using the linear potential. 



30 



1. The p-wave for the bound-free transition 

For the bound-free transition, we can read off from the momentum distributions in Section 
IIV Al that. without the ground state p-wave, 

1 2 

n^ip) = ns{p) = T^V^lni (p) and = n_{p) = . (43) 

Thus, in the scahng hmit, we can introduce the response 

1 2 

Rnop = ^ dpp/;^l ip) (44) 
and then write all the bound-free responses in the scaling limit as 

Rl = -Rnop, Rt = —Rt' = Rnop, RtL' = V^Rnop ■ (45) 

From these equations, one sees that in the scaling limit, the peak position is now identical for 
all four responses, and that the responses differ only by a simple numerical factor. This leads 
automatically to = 1 for all values of q and y. Just like the responses, the spin structure 
functions gi and g2 also peak at y = GeV without the p-wave in the ground state. Apart 
from slight shifts in the peak position, there are no major changes in the asymptotic forms 
of the observables, see Fig. HH For Rl, omission of the p-wave leads to a slight reduction 
in peak height. Rt also has a reduced peak height (a reduction of about 10%). In contrast, 
Rt' maintains its peak height. The only response not to experience any noticeable change 
is Rtl', which was centered around y = GeV to start with, and does not change its peak 
height, either. 

The omission of the p-wave in the ground state simplifies the analytic expressions obtained 
for the responses, and one expects that scaling should set in more quickly. This is observed 
at low q values for Rl and Rt- For Rt' and Rtl', the presence of the p-wave contribution 
plays a negligible role in the scaling behavior. This behavior is consistent with the fact 
that the asymptotic forms of Rl and Rt are more strongly affected by the omission of the 
p-wave. 

For Qi, the results change significantly for low q, see the top panel of Fig. ^1 At g = 
2 GeV, the spin structure function is changing sign at ?/ ~ 0.1 GeV and becomes positive. 
For higher q values, gi peaks at slightly higher y values without p-wave, and the peak 
height is somewhat higher, too. The onset of scaling is very slow, independent of the p-wave 
contribution. 
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FIG. 14: The asymptotic behavior for q — > oo and fixed y of the responses for the bound- free 
transition, without the ground-state p-wave (top panel). We show i?^ (sohd hne), Rt (dotted 
line), Rt' (dashed line), and Rtu (dash-dotted line). The asymptotic behavior for g — > oo and 
fixed y of the spin structure functions g\ and for the bound-free transition, without the ground- 
state p-wave (bottom panel). We show g\ (sohd line) and g^ (dotted line). Note that Ai vanishes 
in this limit, and that A\ = 1. 
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FIG. 15: The spin structure function gi (top panel) and the polarization asymmetry A2 (bot- 
tom panel) are shown without the ground-state p-wave for several values of the three-momentum 
transfer q at fixed y. 

The onset of scaling is dramatically accelerated for Ai by dropping the p-wave, as Rt and 
Rip coincide apart from an overall sign, and Ai is the ratio of these two responses. Therefore, 
Ai scales directly to 1 for even the lowest value of q, if there is no p-wave present. 

While A2 is still scaling slowly, it has quite a different shape now, see the bottom panel 
of Fig. El Instead of starting at its maximum value at large negative y, and decreasing to 
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a minimum around y ~ 0.4 GeV , without the p-wave, it starts out at its minimum value at 
large negative y and steadily increases. 

2. The p-wave for the hound-hound transition 

The situation for the bound-bound transition is similar: for the longitudinal response, we 
observe small changes in peak position and peak height due to the omission of the p-wave, 
just like for the bound-free transition. For the transverse response, we see a more pronounced 
reduction in peak height, as well as the already familiar shift in the peak position. The 
transverse-longitudinal primed response remains largely unaffected by the omission of the p- 
wave, while the transverse primed response has a shifted peak position, and no reduction in 
peak height. These observations are very similar to the observations made above concerning 
the effects of the p-wave on the asymptotic shapes of the responses for the bound-free 
transition. 

The polarization asymmetry Ai again takes the value of 1 immediately, as Rt and Rt' 
only differ by a sign. A2 without the ground state p-wave has the same shape for bound- 
bound and bound-free transitions, markedly different from the shape including the p-wave. 
The spin structure function gi behaves similarly for bound-bound and bound-free transitions, 
too. We do not present any figures for the bound-bound results, as they are so similar to 
the bound-free transition results shown above. 

Strictly speaking, one would have to renormalize the remaining wave function when 
switching off the p-wave. However, the effect of this rescaling will be small, and we omit it 
here. The goal of this discussion was to learn where the p-wave is important, and where not. 
Summarizing, the p-wave has the biggest impact on the polarization asymmetries, where 
switching off the p-wave leads to immediate scaling for Ai and a different shape for A2, 
before A2 reaches its asymptotic value of zero. The spin structure function gi also shows 
sensitivity to the p-wave contribution at low q. 



Our results for the role of the p-wave in Ai confirm the results in j22]. There, the authors 
noted a distinct suppression of Ai at low Nachtmann ^ due to the p-waves. The region of 
large positive where we observe a fall-off of Ai when it is calculated with the complete 
wave function, and find a value of 1 when Ai is calculated without the p-wave contribution, 
corresponds to the low ^ region. 
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E. Model dependence 

In this section, we discuss the influence of the chosen potential on the results for the 
observables. For brevity, we restrict ourselves to the bound-bound transition. Previously 
js^, we found that there is no qualitative difference between the linear, static Coulomb, and 
running Coulomb potential results in the longitudinal and transverse responses. Peak height 
and position are different, but the approach to scaling and the low q oscillatory behavior are 
qualitatively the same. 

First, we discuss the results for the Coulomb potential in the bound-bound transition. The 
responses agree qualitatively with the results for the linear potential, with similar shapes, 
peak heights, peak locations, and scaling behavior. The same holds for the spin structure 
functions gi and g2- For the polarization asymmetry Ai, the decrease from values close to 1 
is slightly more pronounced for positive y values, but the differences are still small. 

This somewhat faster fall-off of Ai is seen more distinctly for the running potential, see 
Fig. El However, the responses and spin structure functions look very similar for the running 
Coulomb potential and the static Coulomb potential. A2 has a slightly different shape for 
the running Coulomb potential, it starts out more flatly at negative y values than for the 
linear potential and the static Coulomb potential. The onset of convergence, however, does 
not seem to be affected noticeably by the potential in the low q region. 

These results for the dependence of the spin observables on the employed model potential 
are consistent with the results found for the behavior of the longitudinal and transverse 
response functions: while the results differ quantitatively, there is no qualitative difference 
due to the different potentials, neither in the results for the highest q valuable attainable 
for bound-bound transitions in our model, nor in the convergence of the results. 

The polarization asymmetries, as ratios of responses, are naturally more sensitive than 
the other observables, but even they only show a minor influence of the potential. 

V. SUMMARY AND OUTLOOK 

We have presented results for responses, spin structure functions and polarization asym- 
metries, calculated within a simple model consisting of a light quark bound to an inflnitely 
heavy anti-quark (or a diquark) without charge. We have investigated the onset of scaling 
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FIG. 16: The polarization asymmetries A\ (top panel) and (bottom panel) are plotted versus 
y for several values of the three-momentum transfer q. The results shown have been calculated for 
the bound-bound transition, using the running Coulomb potential. The results for g = 10 GeV for 
the linear potential and the static Coulomb potential are also shown for comparison. 



and the scaling functions themselves in the bound-free transition, where we are not ham- 
pered by the numerical difficulties involved in calculating the bound-bound transition for 
energies larger than 12 GeV . For the bound-bound transition, we looked at the approach 
to scaling and the duality at low momentum transfer q. 
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We have found that the two polarization responses, Rt' and Rtl', scale roughly like their 
unpolarized counterparts, the longitudinal and transverse responses Rl and Rt- The only 
interference response, Rtl', scales a bit faster than the other responses, but the differences in 
scaling behavior are small. Duality holds qualitatively for the responses, and for the polar- 
ization asymmetry Ai. In fact, Ai scales rapidly, and judging from our model calculations, 
is the most promising observable for the application of duality to extract DIS information 
from resonance region data. These results are reminiscent of the observation that the A^ 
and A^ data show little dependence [g^. 

In contrast, scaling sets in very slowly for A2 and the spin structure functions gi and 
g2. This behavior is independent of the type of scaling variable chosen: convergence for an 
x-type scaling variable and large is as slow as for the y-scaling variable and increasing 
q. We have traced the reason for the slow onset of scaling back to the kinematic factors 
multiplying the responses in the expressions for the spin observables. 

Of course, these are results calculated within a simple model, and the model should not 
be viewed as an attempt at a quantitative description of electron scattering off a proton. In 
nature, it may turn out that the scaling behavior of the responses is different, and might 
compensate in some way for the slow scaling of the kinematic factors. Still, our model 
calculations might provide some useful guidance. It would be interesting to investigate if 
the responses can be separated from the data available at present, and to see if they indeed 
scale and exhibit duality as predicted in our model. 



As in [30'], we have performed calculations for three different potentials. The quantitative 
differences are small for almost all observables, and qualitatively, a different potential does 
not seem to change anything. The one exception is the polarization asymmetry Ai, which 
shows a sensitivity to the employed potential at high, positive values of y. 

We have also studied the influence of the p-wave in the ground state wave function on the 
observables. While most observables undergo only small shifts in peak position and small 
changes in peak height, our results do confirm the larger influence of the p-wave contribution 
on A, at large positive y values observed in Q- However, the other observables seem to be 
fairly robust with respect to the p-wave. 
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